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TO APPEAR IN THE PHYSICAL REVIEW D 

We examine, from both a qualitative and a numerical point of view, the evolution of Kantowski- 
Sachs cosmological models whose source is a mixture of a gas of weakly interacting massive particles 
(WIMP's) and a radiative gas made up of a "tightly coupled" mixture of electrons, baryons and 
photons. Our analysis is valid from the end of nucleosynthesis up to the duration of radiative 
interactions (lO'' K > T > 4 x 10^ K). In this cosmic era annihilation processes are negligible, 
while the WIMP's only interact gravitationally with the radiative gas and the latter behaves as a 
single dissipative fluid that can be studied within a hydrodynamical framework. Applying the full 
transport equations of Extended Irreversible Thermodynamics, coupled with the field and balance 
equations, we obtain a set of governing equations that becomes an autonomous system of ordinary 
differential equations once the shear viscosity relaxation time, Trei, is specified. Assuming that Trei is 
proportional to the Hubble time, the qualitative analysis indicates that models begin in the radiation 
dominated epoch close to an isotropic equilibrium point (saddle). We show how the form of Troi 
governs the relaxation timescale of the models towards an equilibrium photon entropy, leading to 
"near-Eckart" and transient regimes associated with "abrupt" and "smooth" relaxation processes, 
respectively. Assuming the WIMP particle to be a super-symmetric neutralino with mass m„ ~ 100 
GeV, the numerical analysis reveals that a physically plausible evolution, compatible with a stable 
equilibrium state and with observed bounds on CMB anisotropies and neutralino abundance, is only 
possible for models characterized by initial conditions associated with nearly zero spatial curvature 
and total initial energy density very close to unity. An expression for the relaxation time, complying 
with physical requirements, is obtained in terms of the dynamical equations. It is also shown that 
the "truncated" transport equation does not give rise to acceptable physics. 

PACS numbers: 04.40.Nr, 05.70.Ln 



I. INTRODUCTION 



The radiative era of cosmic evolution extends from the 
end of cosmological nucleosynthesis to the decoupling of 
baryonic matter and radiation, covering the temperature 
range 4 x 10"^ — 10^ K (roughly between 1 eV and 1 keV 
During this period cosmic matter can be described |jl|, ||. 
as a mixture of two main non-interacting components 
one, a non-relativistic and coUisionless gas of WIMP's 
(cold dark matter, CDM), the other a tightly coupled 
mixture of non-relativistic baryons, electrons and ultra- 
relativistic matter ("radiation", i. e., photons and neu- 
trinos). The standard approach to the radiative era con- 
sists of using a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) space-time background ||^, ^, ^ whose sources 
are described either by equilibrium kinetic theory 
gauge invariant perturbations [^, or by hydrodynami- 
cal models 0, |8| which in general, fail to incorporate a 
physically plausible description of the matter-radiation 
interaction since they assume a full thermodynamical 
equilibrium throughout the evolution. Since the tight 
coupling between electrons, baryonic matter and radi- 
ation follows from various processes of radiative inter- 
action H, ^ In], mostly involving photons and elec- 
trons, we can ignore the non-interacting neutrinos and as- 



sume the baryon-electron-photon mixture to evolve with 
a common temperature (local thermal equilibrium) and 
to behave as a single fluid, the "radiative fluid". This 
fluid must be dissipative in order to provide an adequate 
macroscopic model for these interactions Ideally, all 
dissipative fluxes (heat flux, bulk and shear viscosities) 
should be taken into consideration in the study of this 
tight coupling. However, in order to deal with a mathe- 
matically tractable problem, while still aiming at a physi- 
cally interesting generalization of previous work, we shall 
study the case in which only shear viscosity is present. 
Bulk viscosity is not significant in the temperature ranges 
we are considering and, although neglecting the con- 
tribution of heat flux carries physical limitations, this is 
compensated by the ensuing mathematical simpliflcation 
of the fleld and transport equations. This approach has 
already been tested in various known and new exact so- 
lutions (12). 

The simplest class of metrics allowing for anisotropic 
shear viscous stresses are the Kantowski-Sachs cosmolo- 
characterized by a 4-dimensional isomet- 
As the source of space-time we consider a 
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gies 

ric group. 

momentum- energy tensor made of CDM (the WIMP gas) 
and the dissipative radiative fluid whose anisotropic pres- 
sure can be identified with shear viscous pressure. Con- 
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sidering the WIMP particle to be the supersymmetric 
lightest neutralino, with mass 100 GeV |l^, |l^, we 
can safely assume that throughout the radiative era these 
WIMP's are non-relativistic, coUisionless and only in- 
teract gravitationally with radiation and baryonic mat- 
ter. It is also reasonable to assume, for the prevail- 
ing temperatures of this era, that the pressure of the 
WIMP gas and the internal energy of the baryons and 
electrons are negligible in comparison with the radia- 
tion equilibrium pressure. Although the radiative era 
is dominated by radiation, the rest mass-energy density 
of the WIMP's is not negligible and dominates that of 
the baryons and electrons, hence the full source, CDM 
plus radiative fluid, can be well approximated by a mo- 
mentum energy tensor in which CDM provides the bulk 
of the rest mass energy ( "matter" , which ends up dom- 
inating the whole dynamics), while the photons ("radi- 
ation") provide the bulk of thermal and dissipative ef- 
fects. The shear viscosity associated with this source 
must satisfy appropriate constitutive and transport equa- 
tions from irreversible thermodynamics that comply with 
causality and stability ||l^, |l9|, |2^, these thermody- 
namical theories are known generically as Extended Ir- 
reversible Thermodynamics (EIT) 
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. The ap- 
plication of such theories to particular physical systems 
requires phenomenological coefficients, like the coefficient 
of shear viscosity, to be provided by Kinetic Theory. In 
particular, for the tight coupling of electrons, baryonic 
matter and radiation and its associated photon-electron 
interaction, the coefficients corresponding to the "radia- 
tive gas" model ^ ^ |27| should be employed. 
The entropy production must be positive definite and the 
relaxation time of shear viscosity must be a positive and 
monotonously increasing function, somehow related to 
the collision times of the radiative processes associated 
with the radiative era. All these timescales must over- 
take the Hubble expansion time as baryonic matter and 
radiation decouple. 

The paper is organized as follows: Sections II to IV 
present and discuss the field equations of Kantowski- 
Sachs geometry for a mixture of CDM and a radiative 
fluid, the application of Extended Irreversible Thermo- 
dynamics and the appropriate set of equations of state 
for the models, as well as the evolution equations for the 
geometric and state variables. The dynamical analysis 
is carried on in Section V by defining a set of normal- 
ized variables, which then leads to a self-consistent and 
well-behaved autonomous system of ordinary differential 
(governing evolution) equations. From the qualitative 
analysis, we identify a saddle point associated with a ra- 
diation dominated FLRW cosmology and contained in 
the invariant set x = 0, associated with the flat Bianchi 
I model. We argue that initial conditions must be de- 
fined near this point. In Section VI we discuss various 
assumptions on the form of the relaxation time for shear 
viscosity, these assumptions lead to the identification of 
a "near-Eckart" and transient regimes, respectively as- 
sociated with a swift and slow rate of transiency. The 



effects of using a "truncated" transport equation are dis- 
cussed qualitatively in Section VII, while Section VIII 
deals with the numerical analysis of the models bearing 
in mind the qualitative results obtained in previous sec- 
tions. The main result that follows from the qualitative 
and numerical analysis is the fact that a physically plau- 
sible evolution is possible only for: (a) initial conditions 
and evolution close to x = 0; ^^^id (b) using the full (not 
truncated) shear viscosity transport equation of EIT. A 
detailed discussion and summary of these results is pro- 
vided in Section 



IX 



II. 



KANTOWSKI-SACHS COSMOLOGIES 
WITH ANISOTROPIC STRESSES 



The simplest non-FLRW cosmological metric allowing 
for anisotropic pressure is that of Kantowski-Sachs (KS) 
models: 



"^dt^+A^{t)dr^ + B^{t) [de^ +sm^{0)d(l)^] . (1) 
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For a co-moving 4- velocity w° , the expansion scalar and 
shear tensor associated with the metric (|^) are: 



(2) 



where a dot denotes the derivative with respect to proper 
time of fundamental observers, which for the KS metric 
(|l|) in co-moving coordinates is given by t, (i. e., A ~ 
A^t = it°j4_a). We consider as the source of (0), the 
following stress-energy tensor: 

T"^ = pM°M^ + vh""" + (3) 

where ft."'' = c~'^u°'v^ + g'^'' and 11" b is the anisotropic 
pressure tensor satisfying HabU^ = H" ^ = 0. The most 
general form of this tensor for the metric (0) is: 



n" fc = diag [ 0, -2P, P, P ] 



(4) 



where P = P{t) is an arbitrary function to be deter- 
mined by the field equations and subjected to an evolu- 
tion law for a given physical model associated with (|l|) 
and (0). The field equations then become: 
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where k = SnG/c^, while the energy balance is given 
by: 



p + {p + p)e + 6(7P = o. 



(8) 



A. Mixture of cold dark matter and a radiative 
fluid 

We will assume that the stress-energy tensor (|^) corre- 
sponds to a mixture of a non-relativistic gas of WIMP's 
and a radiative fluid with shear viscosity corresponding 
to a "tightly coupled" mixture of photons, electrons and 
baryons sharing a common temperature T. Hence, p 
and p in eqn. (|^) are the total mass-energy density and 
equilibrium pressure given by: 



P = Pw + /Ob + Po + P-,, 

p = p, + Pb + p. + P-,, 



(9) 



with the subindices 7, w, 6, and e denoting photons, 
WIMP's, baryons and electrons, respectively. The three 
latter components satisfy each the equation of state of a 
non-relativistic ideal gas: 

p„ = (m^.c^ + 3fc^T'„/2) n„, p„ ~ k^T^, 
Pb = {m^c^ + 3k^T/2) n^,, p^, = k^T, (10) 
= (m^c^ + 3k^T/2) n,, p, k^T, 

where m„, rrib, are the respective particle masses 
of the WIMP's, the baryons (a proton mass) and the 
electrons, fc^ is Boltzmann's constant, T„ is the tem- 
perature of the WIMP gas and T is the common tem- 
perature of the radiative mixture. During the radiative 
era creation/annihilation processes cease to be significant 
and so the particle number densities, n„, rib, n^, satisfy 
conservation laws of the form: 

•h + nO = 0, with n = n^., n^,, (11) 

which can be integrated, leading to 

N 



with 



N 



N„, iVb, (12) 



where iVb, are the constant {i. e., conserved) 

number of WIMP's, baryons and electrons, respectively. 

The radiation component of the radiative fluid can be 
given either in terms of: 

(«) the Stefan-Boltzmann law: 

Pf = a T\ = a TVs, (13) 

where a = tt^/c^ /(ISfi.'^c'^) is Stefan-Boltzmann constant, 
or 

(m) an ultra-relativistic ideal gas: 
p- = 3 fc^ T, = k^ T, (14) 



where is the number density of ultra-relativistic par- 
ticles, subjected to a balance law analogous to (^ and 
given by an expression similar to ([ij) with the conserved 
photon number N^. 

In order to simplify eqn. (^), we can examine the ra- 
tios of particle numbers and rest mass densities of the 
different particle components. Considering the currently 
estimated H ratio of photons per baryon, we have: 



2.67 X IQ-^n^ h^, 



(15) 



where f2b is the baryon abundance today (roughly 
0.04 ± 0.01) and /i ~ 0.7 is the adimensional Hubble 
factor ||2^. Regarding the WIMP gas, if we assume that 
it is made up of the lightest supersymmetric neutralinos 
with m„ - 100 GeV |l|, |l6[ 0' have that |l|: 



2.82 X 10"^ h^. 



(16) 



where Q,^ ~ 0.3 ± 0.1 is the neutralino (CDM) abun- 
dance today. Using eqns. (p^, (|l2|) and (p^), we can 
rewrite (0) as: 



m„c n„ \ \ + 



Vb , , v^T^ 
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p = n^k^T ( 1 + i^b + 
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Bearing in mind that for electrons N-^,, so that 

j^b, while 4; nib ~ 1 GeV, ^ 1, and 
j/e ~ t^b ^ 1, then. 



10-'^ « 1, 



while for the temperature range 4 x 10"^ K < T < 10^ 
K, we have: 



0.013 < ' < 3.25, 

Tl-y rZ ^ J. 



(18) 



showing that the radiative era is initially radiation dom- 
inated but rest mass energy density is not negligible and 
ends up becoming dominant. Therefore, even if T^/T 
in eqn. ( |l7|) is not negligible, we have: 



m„c n„ -|- 3n^k„T, 



p ~ n^k^T, (19) 



The same type of approximation can be obtained if we use 
the Stefan-Boltzmann law ([l3|), since the ratio of pres- 
sures in eqns. ( p^ and (p3|): pt,''/Pb = ciT'^ /[3ni,kg], is 
proportional to z^b (likewise for WIMP's). Therefore, 
the mixture of a gas of WIMP's and a radiative fluid can 
be accurately described in the desired temperature range 
by the approximated equation of state: 

p = mc^n^") + p^^'l p = p''') = p(''V3 
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with 



= p„ (20) 

where p'''-' follows from either one of eqns. ( p^ ) or (14), 
hence n^"^"^ — and p*-™-* = m^c^n^. For the remain- 
ing of this paper, the superindices and will 
refer to quantities associated with photons ( "radiation" ) 
and WIMP's ("matter"), respectively. From eqns. (|^) 
and (^), the equation of state ( pO| ) can be given as the 
following constraint: 
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0. (21) 



III. 



EXTENDED IRREVERSIBLE 
THERMODYNAMICS 



If the source (|^) is meant to describe a mixture of 
non-rclativistic CDM and a radiative fluid (as argued 
in previous sections), the anisotropic pressure must be 
identified with a shear viscous stress of the latter fluid 
and must be compatible with a suitable thermodynami- 
cal formalism. We shall consider the so called "Extended 
Irreversible Thermodynamics" (EIT) (2^, |2l|, |2|, |2|, a 
theory complying with causality and stability require- 
ments ||2^ and supported by the Kinetic Theory of gases, 
Information Theory and by the Theory of Hydrodynam- 
ical Fluctuations When shear viscosity is the only 

dissipative agent, the corresponding generalized entropy 
current, S"", obeying the usual balance law with non- 
negative divergence, and up to second order in II"'', 
takes the form: 



nSu", 



2nT 



(22) 



where 5* is the entropy per particle, n 
{l + i^„) n^^^ is the total particle number density (« n^-^^), 
a is a phenomenological coefficient to be specified later 
and S'(e) is defined by the equilibrium Gibbs equation: 



nTS, 



P- {P + P^''')- = 3p + 4pe = - 
n 

where we have used eqns. (O) and (E3) to eliminate p. 
Fulfillment of the second law of thermodynamics requires 
5'° -a > 0, which from the definition of 5° and S in 
( p^ ) leads to: 



together with 



S > 0, 



2 77' 



(24) 



(25) 



and the evolution equation of the viscous pressure, i. e., 
the transport equation |e3, E9|: 
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(26) 



where 77, T, t„i, aab, ^ab are the coefficient of shear 
viscosity, the temperature, the relaxation time, and the 
shear and shear viscosity tensors, respectively. The pa- 
rameter eo can take only two values: eg = 1 ("full" 
transport equation) and eo = ("truncated" trans- 
port equation, also known as the Israel-Stewart equa- 
tion |l|, ^), while Eckart's non-causal transport 
equation follows by setting r„i — 0. The coefficient of 
shear viscosity as well as other related quantities can be 
obtained by a variety of means ^ including Kinetic The- 
ory, Statistical Mechanics or both [|5[ ^ |3^ . Unless 
specifically stated otherwise, we shall consider only the 
full transport equation cq = 1. We will discuss the im- 
plications of the truncated equation (eo = 0) in Section 
VII. Evaluating S from (p^ and using eqns. ( p3| ^ 
and (I2I) we obtain: 



S 



2 r; n T 



2 [ 5(e) - S 



> 0, 



(27) 



where we have used eqn. (p^ ) and eo = 1 in cqn. (|2q). 
Notice how S takes a very simple form, illustrating the 
role of r„i as the reference timescale associated with the 



,) > S. Also, the second 
ll)) is fulfilled if ?7 > 
> S hold. Another 



entropy change from S to S'( 
law of thermodynamics, (i. e., 
or, equivalently, if t,^i > and ^(e) 
important requirement that follows from the second law 
of thermodynamics and the stability of equilibrium states 
is that 5 be a convex functional, i. e., S < 0. For 
the KS models all quantities depend only on time and 
so a necessary (but not sufficient) condition is given by 
5 < 0, which leads to: 



1.^ 



Uab n° 



2 77 (7ab Tl"" > 0. 



(28) 



For the applicability of the general relations, eqns. ([2^), 
(H), (H), (H), (|^ and (H) to the models we are 
concerned, we must impose the equation of state (2C), 
with either one of the choices eqn. (|l^) or eqn. (14). 
As a physical reference to infer the form that the co- 
efficients T„,, 77, a may take for the radiative fluid, 
consider the "radiative gas" model associated with the 
photon-electron interaction |2^, ^ , characterized 

by p, p complying with the equations of state discussed 
in section III. For the radiative gas the forms of 77, a 
in terms of the relaxation time of the dissipative process, 
T„i , are: 
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(29) 



where p^^'^ is either or n^^^k^T (eqns. (|T^) or 

(p^), respectively) and the subscript ''{rgy emphasizes 
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that these quantities are specific to the radiative gas. 
Applying eqn. into eqns. ( p2|) and (p3|), we get for 
the entropy per particle: 
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p(r) nW T' 



3 CT P 



(30) 



(31) 



where we have neglected the entropy of the non-relativ- 
istic matter, so that n^™) + n^'') = (1 + v„) n^'"> « 
and S is approximately the photon entropy (this is jus- 
tified because the number of photons is so much larger 
than that of WIMP's, baryons and electrons). The trans- 
port equation (pi) becomes: 



P 



° (r) 
— p^ '(7 

5^ 
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p 



Act 



P 



(r) 



p2 



0, (32) 



where A = 1 + l/(2Ao) and Aq = 1/2 or 2 (for 
given by the Stcfan-Boltzmann law (O) or the ideal gas 
law (H), respectively); also condition (|27D takes the form: 



0). However, since both EIT and Eckart's theory assume 
near equilibrium states, quantities like P^ and aP, , 
appearing in eqns. (^) to (^ must be small, hence 
the ratio T'^/n'^'"^ is nearly constant and so we can also 
assume that S'(e) given by eqn. (34) is approximately 
valid for the ideal gas law. 

Conditions ( ^ ) and 6^ S < associated with eqn. 
( p8| ) must be satisfied by any self consistent thermody- 
namical system. The importance of these conditions will 
become evident when discussing the numerical integra- 
tion of the evolution equations. The relaxation time, 
T,„,, is qualitatively analogous to and larger than the 
mean collision time between particles and it may, in prin- 
ciple, be estimated by collision integrals provided the in- 
teraction potential is known. Since we are concerned with 
mixtures of baryons, electrons and photons evolving in 
the temperature range 4 x 10'^ K < T < 10^ K (from the 
end of cosmic nucleosynthesis to decoupling) , convenient 
references for comparing t,„, are the collision times as- 
sociated with Compton and Thompson scatterings [^ : 
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^7 J 



(36) 



S = 



15 P2 



2 



S 



4 pi"-) T 
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Notice that S < S'(e) and t„i > must hold in order for 
5 > to be satisfied, while (|23| ) implies that we must 
have a P < in order that S(e) > 0. Regarding the 



interpretation of 5(e) , if we assume SVp\ > 



and 



the Steffan-Boltzmann law dl^), eqn. (23) then yields: 

4 a 



3 n(-) ' 



(34) 



a function that is only constant (equilibrium photon en- 
tropy) if P = 0. This constant is given explicitly by the 
black body formulae psj: 



30C(3)aT3 



45C(3) 



3.60 fc. 



(35) 



where is the Riemann zeta function. However, if we 
characterize the evolution to equilibrium as P — > 0, 
then, as this evolution proceeds, S(^e) in eqn. ( p^ must 
tend to the constant entropy given in (|3|). Hence, we 
can identify eqn. ( p5[ ) as the equilibrium state associated 
with eqn. (p7|) and eqn. (B3h, attained as P ^ 



and S — > S'(e) {i. e. as the radiation relaxes) in the 
timescale provided by the relaxation time t„,. For the 



ideal gas law, equation (g3|) does not yield eqn. ( |34D , but 



5, 



ln(r^/n'^''^), an expression that coincides 



with eqns. (|3J) and ( pS] ) only in equilibrium (if S(^e 
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1/2- 



3/2 



, (37) 



where cr^ , Bq, m„ , and h are the Thompson scatter- 
ing cross section (6.65 x 10~^^ cm^), the hydrogen atom 
binding energy (13.6 eV), the electron mass and Planck's 
constant, respectively. Equation (|3^) is obtained from 
the number density of free electrons provided by the Saha 
equation. Notice that we are using the baryon number 



density, 



and not the number density of WIMPS, 



^ higher temperatures in the range of interest, 
Compton scattering is the most efficient radiative process 
keeping baryonic matter and radiation tightly coupled, 
though it is no longer effective in lower and intermediate 
temperature ranges (T < 10^ K). The photon-electron 
interaction of the radiative era requires that microscopic 
collision times i^, tc-, as well as r„i, be much smaller 
than the timescale of cosmic expansion given by the Hub- 
ble time, approximately = 3/0. For the lower 
temperature range of the radiative era, just before re- 
combination, Thomson scattering becomes the dominant 
radiative process, so that the decoupling of baryonic mat- 
ter and radiation can be associated with the condition 
= t^, which should be approximately equivalent to 



IV. EVOLUTION EQUATIONS 

Since we need to determine a self-consistent set of or- 
dinary differential equations governing the KS models, it 



6 



is convenient to express the field equations and eqn. (|2l| ) 
in terms of Q and a by eliminating A, B from (0): 



A e 



— = — + <T 



which leads to: 
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kP = (T 



while, using eqns. (|3 
( pT| ) becomes: 
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a& 
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140), the equation of state 
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(38) 

(39) 
(40) 
(41) 
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Kmc' n*^™^ 



0. (42) 



We can eliminate B from the equations above with the 
help of equation ( p9[ ) (which becomes a constraint) and 
the second of equations (|38| ) (the evolution equation for 
B). Using equation (pO[), equations ( p] ) and ( p2| ) then 
become: 

a + a' + aQ - ^ + '^mc'n^"'^ + np^''^ - kP = 0, 

(43) 
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3 K p' 
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0, (44) 



which are the evolution equations for cr and O. Since we 
are assuming particle number conservation, an evolution 
equation for non-rclativistic particle number density (the 
WIMP's) follows from equation (11): 



+ e ^ 0, 



(45) 



a conservation law satisfied also by n'^''^ = n^™'/i/„ (if 
using the ideal gas law ([T^)). Another evolution equation 
is provided by (|g), which applied to the equation of state 
(p^ and using <Jab = 6 cr F yields: 



•(r) 



+ o P^''' 6 + 2 cr P = 0, (46) 
o 

becoming the evolution equation for p'^''-' . The transport 
equation ([33) derived in the previous section, namely: 
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8 (r) 
-p^ >C 
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4^P2 ^ (47) 



is the evolution equation for the shear stress P. In addi- 
tion to these evolution laws, eqn. (33) can be thought of 
as an evolution equation for S, while we can transform 
eqn. ( [46|) into an evolution equation for T by using 
either ^13| ) or (p^), leading to: 



T 
f 
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3" 



cr P 



Ao P 



(r) 



= 0. 



(48) 



Equations (|43|), (||), (||), (g), and (|3) represent a 
self-consistent and closed system of first order ODE's for 
7T,(™), P, p^'^\ cr, 0. Notice that this set of evolution 
equations is fully determined if r,„i is known. In the 
cosmological context r,„i might be proportional to the 
timescale defined by the expansion scalar (approximately 
the Hubble time): 



r,„i oc t 



3/e. 



(49) 



Alternatively, and depending on the temperature and en- 
ergy range one is considering, r^ei could be identified as 
proportional to a collision time (say, Thomson or Comp- 
ton scattering) given by equations (pm or (p7|), i. e., 



r,„i oc (n, T) , 



a tc(n,T), (50) 



Since we can assume in equation (|4^) two different equa- 
tions of state for the radiation component, strictly speak- 
ing, equations (^)-(p7|) constitute two different systems 
of evolution equations parameterized by the two possible 
values of Xq and A: 

-Stefan-Boltzmann law: 

3pW = aT*, with Ao = 2 and A = 5/4, 
-Ideal gas law: 

= n^^'h^T, with Ao = l/2 and A = 2. (51) 

V. DYNAMICAL ANALYSIS 
A. The governing equations. 



Let us define: 



Q 



p 



M' 



(52) 



a ratio that is related (from eqn. (pO[)) to the deviation 
from equilibrium: 
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where: 



An = 



P 



(r) 4 C'(e), 



^ S(p\, Stefan-Boltzmann law 



Ideal Gas law 



(53) 



with S'(e) given by eqn. (p4|). Therefore, Q must be 
a small quantity for the cosmic times we are interested 
in. Note that when P (and hence Q) is zero, the 
shear vanishes and the Kantowski-Sachs model reduces 
to an isotropic FLRW model. In principle, P (and 
hence Q) can be positive or negative, but since P 
represents a viscous pressure it should be negative in the 
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expanding regime. We shall focus henceforth on the case 
in which Q is negative. The energy conditions imply that 
— l<(3<l/2. In addition, on physical grounds we 
expect the second term in equation (jg) to dominate the 
third term in this equation, which is satisfied whenever 
cr^(5^/(8/3)^ < 4. However, since physically we expect 
Q"^ to be small, these constraints are satisfied handily. 
We introduce now the new normalized variables: 



6/3' 



„2 (m) 



(m) 



62/3 



(r) 



up 



(r) 



(6/3)2' 



(54) 



a definition that is motivated by the fact that 6/3 is 
approximately the Hubble expansion factor H ~ 6/3 + 
(Jabn°"n'^ for unit space-like vectors defined by n^Tij, — 
1, u°'na = 0, hence r2(m) and f2(r) are approximately 
equivalent to the observational parameters (the fi's) for 
the CDM and radiation contents of the mixture. We also 
define (for 6 > 0) the new independent variable: 



Eventually, the models re-collapse and 6 changes 
sign. At the point of maximum expansion (when 6 = 0) 
the variables above diverge and the normalized equations 
are no longer valid. However, for the times we are inter- 
ested in, in the expanding phase far from re-collapse, the 
above variables and equations are valid. Indeed, in prin- 
ciple we can use the above system to follow the evolution 
of the models all the way back to the big bang. From 
equation (H) we can see that in this regime the curva- 



ture is small and that the variables E , i7 



(m) 



(r) 



are 

well-behaved. Compact variables can be defined by nor- 
malizing with 6^ + B~'^ (instead of 6^) that are valid 
for all times ^2\; however, the physical assumptions used 
here are not valid at later times. 

In addition to equations ( p6| ) to (^9|), evolution equa- 
tions for 71, T, S(^) and S follow by using the variables 
defined in eqns. ( |5^ ) and ( [55| ) in equations (45), (|4^), 
dH) and (H): 



7i(t) 



= —3 n, 
where: 



(62) 



d 
di 



6 d_ 
"3 d? 



, 

3 *' 



(55) 



the evolution equations ( ^3| ) to ( ^ ) become: 

= - n^m) [ 1 - 4s2 - r!(„,) - 2n^r) ] , (56) 



T' 



-T 



S Q 



Ao 



'S'(e) — — 3 Ao S Q, 



(63) 



(64) 



- E 

5 



Q 



r,., 6/3 



(A - 2) E Q2 



(57) 



^[r) = ^^(r) [2 + 2E (Q — 2E) — rJ(m) ^ 2f2(r)] , 

(58) 

E' = (1 - 2E) (1 - E2) - n(,„) (2 - E) /2 - 

(1 - g - E) , (59) 



where we have used 

36 _ 6^ _ 
"62" ~ e" " " 



1 - 2 E2 



(60) 



and a prime denotes differentiation with respect to r. 
Note that this last equation implies that 6 is monoton- 
ically decreasing. The constraint ( ^2| ) becomes: 



(m) 



(r) 



X = 1 - E2 

We note that: 

X' = X [ 2E (2E - 1) + 



52 e2 



(61) 



2n 



(r) 



Clearly x — is an invariant set of the above differen- 
tial equations, which corresponds to the Bianchi I (zero 
curvature) sub-case. 



S' 



2t, 



5(e) - ] = 



15 t 



4 T„, 



^ Ao 



(65) 



where Ao is given by eqn. (|5^), the subindex i denotes 
evaluation at an initial time r = t^, and — 3/6, is 
approximately the Hubble time which follows from eqn. 



= Q-exp 



1 + 2E 



2 , ^(™) 



(r) 



dr}, (66) 



while the relation between physical time t and r follows 
from and (Ml: 



t — I dr , 



(67) 



where is given by (|66|) above. 



B. Qualitative Properties 

Consider the dynamical implications of assuming that 
T„i is given by eqn. (E^), namely: 



7o 
6/3 



70 t,_ 



(68) 
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where 70 > is a constant. For the range we are 
interested in, 70 < 1- Equation becomes: 



= - (A-2)I]Q^ 

5 70 



(69) 



and so equations (|5^), (58), ( p9[ ) and ( |6^ ) now consti- 
tute a closed four-dimensional system of first order au- 
tonomous differential equations for (Jl(r'), ri(„i), S, Q), 
a system that depends on the value of the constant pa- 
rameter 7o. Moreover, from equation (|6l|) and the above 
discussion, in the regime we are interested in 
and S are bounded and physical conditions imply that 
— 1 < Q < 1/2. Consequently a local analysis of the sta- 
bility of the equilibrium points of this system will provide 
useful dynamical information. 



Setting the right-hand-side of equation (69) to zero we 
obtain: 



(A-2) s g2 



Q 

70 



8S 



- 0, 



(70) 



which is a quadratic equation for Q if S is given and 
A ^ 2. 

Stefan- Boltzman law. (A — 5/4): 

The equilibrium points at finite values are (note that all 
such points are given below; however, (ri(m), S^) are 
not necessarily bounded by unity): 

(i) ri(„i) = = six points: 

(ia) S = 1/2 and 

Q± = 4 [370 ± + 2/15 ]/(97o'); 

(ib) E ^ 1 and 

Q± = 2 [370 ± + 8/15 ]/(97o'); 

(ic) S = — 1 and 



Q± = -2 [870 ± V97o' + 8/15 ]/(97o'). 
(ii) — 0, fi(r) / 0, four points: 

(iia) r!(,) = -(1270 + 17/15)/(972), and 



■8/15/(370); 



=Q± = ±2^3^ 

(uh') O, . , ~ (4-H37o)[2+7o(l±")] 

= (4 + 37n)/[2 + 37o(lTu)], and 
0± = (2 ± 37ou)/(4 + 370), where 



(iii) n 



(r) 



u = 2v/l + (87o/5)(4 + 37o)/(37o); 
J = = Q, fi(r) = 1 is a particular solution. 

0, ^l(rn) 7^ 0, one particular solution: 







Q, and 



1. 



(iv) ri(„i) 7^ 0, f2(r) ^ 0, two points: 
17(„.) = (40 + 79/7o)/(647o), 

r!(,,)--15(8 + 37o)/(1287o), 



±V5(8 + 37o)/(27o)/8, and 
= tV2(8 + 37o)/(457o). 



The two equilibrium points given by (ia), can be shown 
to lead to a value for Q which is unphysical in the 
sense that its magnitude is much too large (in fact, when 
^(m) = = Q would be expected to vanish). The 

equilibrium points (iv) are also unphysical: since clearly 
E^ < 1/4 and we get from (iv) that E^Q^ ^ 1/4^ this 
would imply that > 1 ! At the equilibrium points (ii) 
when fl(^r) 7^ 1 we have that E = 20Q/[7o(15(52 - 32)], 
which leads to a quartic equation for Q; however, for 
physical values of the parameter 70, this equation has no 
real roots and hence no solutions of physical interest. The 
two equilibrium points given by the particular solution of 
(ii) and by (iii), corresponding to FLRW models, namely 
{^{r),^{ni),^,Q) given by (1,0,0,0) and (0,1,0,0), can 
easily be shown to be saddles. The equilibrium points 
(ib,c), namely (ri(r), E, Q) = (0, 0, ±1, (3±), which 
belong to the invariant set x = (i. e., correspond to 
a Bianchi I model with no curvature) have eigenvalues 
3, -1/70 + 3A2EQ/2, 2 - 2Eg, 4 - E. . In this case, 
the eigenvalue — 1/70 + 3A^E(5/2 can only be positive 
if we take the 'positive square root' of eqn. (37), (i. e., 
Q = (3+ ). However, for {A^,7o}<l, i. e., X^^o<l, 
it follows that the eigenvalue 2 — 2E(5_|_ can never be 
positive. Consequently, this equilibrium point can not be 
a source. 

We note that there are no sinks at finite values. How- 
ever, this is to be expected since the models evolve to- 
ward maximum expansion at which the variables become 
unbounded. The models subsequently re-collapse. 
Ideal Gas case (A = 2): 

Equations (|5^), ( p8| ) and (|5^) remain unchanged and 
the governing equation (|69|) becomes: 



5 70 

At an equilibrium point we immediately see that: 

8 70 E 



(71) 



(72) 



However, there are no major qualitative changes in the 
analysis. In particular, there are fewer equilibrium 
points: 

(i) r2(m) = = three points: 

(ia) E = 1/2 and Q = -47o/5; 

(ib) E = 1 and Q = -870/5; 

(ic) E = -1 and Q = 870/5. 

(ii) f7(„i) 0, r2(r) / 0: 

E = = Q, = 1 is a particular solution; 

otherwise, 

A quadratic equation for S which gives an- 
other two equilibrium points. 

(iii) r2(r) = 0, ri(„i) ^ 0, one particular solution: 

E = = Q, and = 1. 
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(iv) ri(„i) 7^ 0, f2(r) 7^ 0, two points that imply: 



A. The "near-Eckart" and the transient regimes. 



I]2 = 5/(1670) > 5/16 (for 70 < 1), which 
leads to fl(^,n)+'^^{r) = l-5/(47o) < 0. 



ASSUMPTIONS ON THE RELAXATION 
TIME 



VI. 



In the general dynamical system, eqns. ( |56| ) - (|59|), the 
relaxation time r,ei needs to be specified in order for the 
governing system of equations to be closed. For the qual- 
itative dynamical analysis we have assumed (|6^) , so that 
1/(^,010) is a constant, since otherwise we would either 
not be dealing with an autonomous system, or would be 
looking at an autonomous but much more difficult dy- 
namical system. However, equation (68) is a simplify- 
ing assumption that cannot be supported by thermody- 
namical arguments, perhaps a more realistic assumption 
would be to consider instead: 



7(t) 

e/3 



= 7(t) t„ 



(73) 



so that equations (|57|) and (|65|) become 
8 E Q 



7(t) 



- (A - 2) E 



S' 



7(r) 



S 



(74) 



(75) 



where 7 = 7(t) is a function that could be suit- 
ably adjusted so that t„i has a form that is qualita- 
tively analogous to that of microscopic timescales like 
eqn. ( ^6|) or (^7|), timescales that are physically rele- 
vant for the matter source under consideration. The ratio 
3/(r„i0) = i„/T„i = 1/7 should provide a comparison 
of the timescale for the relaxation (transient) effects in 
the radiative fluid with the timescale of cosmic expan- 
sion. Hence this ratio should approach unity as baryonic 
matter and radiation decouple, so that 7 w 1 should be 
a consistent choice for near decoupling conditions, while 
7 :§> 1 or 7^1, correspond to after decouphng (late 
times) and much before decoupling (earlier times). Ide- 
ally, we should obtain r,„i from collision integrals as- 
sociated with each of the various radiative processes oc- 
curring in the radiative era, but such an endeavor would 
merit a separate paper by itself and will not be attempted 
here. Instead, we will consider t„i as an "effective" re- 
laxation time, encompassing the different radiative pro- 
cesses. We will examine the non-transient limit (or near- 
Eckart regime) and use the dynamical equations them- 
selves in order to suggest a suitable form for r,„i . The dis- 
cussion of this section will be complemented and tested 



In order to examine the relaxation process as ^ 
S'(e), we will assume that 7(t) in (^3|) is a smooth 
function so that we can always expand it in the form 
7 w 7(0) + 7'(0) r + 7"(0) tV2. Therefore, at early 
times T « we can always associate the constant 70 in 
eqn. ( |6^ ) as 70 = 7(0), so that the corresponding form 
of T,„i is approximately correct at least near t = 0. 
Also, at early times we must have 7(0) < 1, and so 
1/70 ^ Ij while E(0) and Q{0) are necessarily small 
quantities. Hence the terms E and SQ^ in eqn. ( [T^ ) 
will be much smaller than (5(0) /70 and so we have that 
near r = 0, Q' « — Q/70, so that: 



Q « Q(0) exp 



70 



S - ^(e) - ^ Ao Q\0) exp (— 
8 V 70 



numerically in section VIIl (see also figures 5 and 6). 



(76) 



where Aq is given by eqn. (|53|). Since the process 
of relaxation to equilibrium can be characterized as the 
decay of the dissipative flux Q ^ as S grows and 
asymptotically approaches the equilibrium state given by 
eqn. (35), the numerical value of 70 in ( |76| ) may be 
interpreted as a measure of the "rate of transiency" in 
terms of how fast or slow the system accomplishes this 
relaxation in comparison with the timescale provided by 
. We can then identify two possible situations 

• The near-Eckart regime. If 70 <^ 1, then eqn. 
( [76|) indicates a very fast relaxation with a very 
abrupt decay of Q and S to S^e)- In the 
very limit 70 ^ (so that t„i — > as well) 
we have Q — *■ Q(0) S{t), so that the relaxation 
is inflnitely fast, in agreement with the non-causal 
nature of Eckart's classical theory. The fast relax- 
ation associated with a very small 70 implies a 
very short duration of the relaxation process (small 
r,^i), indicating that the approximation 7 = 70, 
as well as the expressions in eqn. ( fz^ ) for Q{t) 
and S{t), are approximately valid for the whole 
evolution time. Hence, since Q practically van- 
ishes very quickly, we have for most of the evolu- 
tion time that S ~ 'S'(e)5 agreeing with the fact 
that in Eckart's theory the entropy is rigorously 
and unambiguously given by its equilibrium form 
(the hypothesis of "local equilibrium", see section 
1.3.1 of ll^). The near-Eckart regime is appropri- 
ate to describe a given radiative process for which 
microscopic timescales are much smaller than the 
cosmological expansion timescale t^. See figures 
5. 

• The transient regime. If 70 < 1 but 70 ~ 
O(lO-i) - 0(1), then Q and S also decay 
but the relaxation process is much slower, hence 
the term "transient" . In this case, the expressions 
in eqn. (|7^) and the assumption 7 = 70 are only 
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good approximations for 7, Q{t) and S{t) near 
T = 0. In general, we must use: 



Q w Q(0) exp 



dr 

7(t) 



S ^ S, 



15 



Ao g^(0) exp 



leading to eqn. 



S" 



2 dr 

7(t) 

51) and to: 

2 (1+7') (5(e) -5) 

7^ 



(77) 



(78) 



so that the fulfillment of S" > and S" < 
can be examined in terms of 7(t) and S'(e). 
Sufficient (but not necessary) conditions follow by 
demanding that 7 is a monotonically increasing 
function (7' > 0) and S'^^-) > 0, S['^^ < 0. Another 
condition on 7 is furnished by eqn. (|28|), which 



together with eqns. 
& yield 




< 2 

7 



5 



7'g' 
1 



21]^ 



f7 



(m) 



r2 



(r) 



(79) 



a condition that should be tested numerically for 
any given choice of t„i. 

The relaxation time r^oi can be approximated by a 
given 7 = 7o, but as the decoupling era is reached, 
7 must increase to 0(1) allowing for t„i to 
overtake t„. 



B. Dynamical relaxation times. 



Let us consider the following ansatz: 



1 

7 



8 E 
'5 Q 



[1 + C(r) 



(80) 



which provides an exact relation for the relaxation time 
in eqn. ([73|). In regimes in which the Eckart theory is 
a good approximation, we can assume 7 « 70 ^ 1, 
so that eqn. (|8^) imphes Q cx S but jS/Q] > 1 
and C is 'small'. Such regimes would correspond to a 
radiative process that takes place in timescales smaller 
than a mean collision time, thus decaying very fast to 
equilibrium. In case we wish to consider processes tak- 
ing place on timescales comparable and larger than main 
collision times, then transient effects are important and 
the near-Eckart regime is no longer appropriate. In par- 
ticular, we can construct an expression for t,„i that 
acts as an 'effective' relaxation time that encompasses 
the relaxation times for the main radiative processes act- 
ing in the radiative era under consideration. Having this 
idea in mind, a reasonable and more general expression 
for C can be obtained from the dynamical equations 



themselves. Consider the conditions discussed above for 
a near-Eckart regime, and assume that: 



= -Mo S, 



(81) 



where fiQ is a non-negative constant. We obtain an 
expression analogous to eqn. ( pO| ) by substituting eqn, 
(31) into the (full EIT) evolution equationi 
(59), leading to the consistency requirement: 



1 

7 

where 



70 



60S' 



1 - MO - ^ ) X, (82) 



1 

70 
ao 



Mo 



5mo 
Mo, bo 



> 1 



1 + (A - 3)mo, (83) 



^-iid and x is the curvature term 



X 



1-n 



(m) 



(r) 



given by eqn. (|61|), a term that is very small and can be 
neglected. This is an extremely simple dynamical rela- 
tion for the relaxation time, and for a wide range of con- 
ditions it might be a very good approximation. In addi- 
tion, it has some important physical properties. For early 
times in which x/S, ^(m)-, S^ are very small, we 
have a near-Eckart regime associated with 7 w 70 <^ 1, 
as expected. This is appropriate for the Compton scat- 
tering, the dominant radiative process in the early part 
of the period under consideration, a process that quickly 
thcrmalizcs and ceases to be effective. At later times, as 
ri(m) increases toward a value of order unity at recom- 
bination (and to a lesser extent, the curvature term also 
grows), the constants 70, ao, 60 in eqn. (p3|), as well 
as the initial conditions, can be selected in such a way 
that 7, , given by (^, increases sufficiently as to allow 
Tr^i to overtake and to approach the characteristic 
timescale of Th omson scattering in eqn. (37). We show 
in section VIII that adequate parameters and initial con- 
ditions can be found so that r,oi associated with eqn. 
(^) has the expected behavior (see figures 6). 



VII. TRUNCATED THEORY 

The discussion so far has been based on the full trans- 
port equation. In order to appreciate the effect of con- 
sidering the truncated transport equation, it is useful to 
rewrite equation ( p6| ) with 77 and r,ei given by eqns. 
( p9[ ) and (^) in terms of Q and S, thus allowing 
the full and truncated equations to appear jointly. This 
yields: 



8 S 



1 



7(r) 



4(1 -eo) Q ~ (Aeo-2)SQ2. 

(84) 

where the full and truncated theories are respectively 
given by eg = 1 and eg = 0. Assuming a transient 
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regime, instead of eqn. (uq) we have near t = 0: 



Q(0) exp 
15 



70 



+ 4(l-eo) T 



(e) 



8 



AoQ'(0)exp<^ 2 



7o 



+ 4(l-eo)T 



(85) 



so that using the truncated transport equation (eq = 0) 
introduces a large linear term (oc 4r) that is absent in the 
full theory. This linear term could change dramatically 
the form of Q and the relaxation of S to 'S'(e). As we 
show below, the truncated transport might lead to Q S 
being positive (implying that ^(e) < 0), so even at this 
level (early times) there might be problems of consistency 
in the truncated approach of EIT. 

A comprehensive analysis in the case of the truncated 
theory, similar to that presented in Section V, can be un- 
dertaken. The evolution equations (^6|), (|5q), and (|59| ) 
remain unchanged, while the evolution equation ( |57| ) for 
Q must be replaced by the truncated equation that fol- 
lows by setting eo = in eqn. (p^) 



Q' 



70 



At an equilibrium point: 



2 S Q2 



1 

70 



- 4 g + 2 E Q2. (86) 



0. (87) 



Q 



However, we immediately note from equation ( p7|) that 
close to an equilibrium point for small S and Q , 



S ~ - 4 
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so that for 70 > | we have that E and Q have the 
same sign (unlike the non-truncated theory case), and 
hence it is immediately clear that there will be a different 
qualitative dynamical behavior in the truncated theory. 

In particular, close to the FLRW equilibrium point 
(ri(r), r2(„), S, Q) = (1, 0, 0, 0) (which we will use to de- 
termine the initial conditions in our numerical analysis 
in the following section), from a calculation of the corre- 
sponding eigenvalues we have that: 



(m] 



Gc e 



S, Q 



(1 — fi(r)) cx 



-/3e° 



(89) 



(90) 



from equations (56) - ( p9[ ) in the non-truncated case, 
where /? is a constant, and 



a± 



1 
2^ 



70 - 1 ± V (70 - 1) - 470 (1 + 870/5) 



(91) 

For physical values 70 < 1, both a± have negative real 
parts and as noted earlier, this FLRW equilibrium point 
is a saddle (e. g., for 70 = 1/2, a± = -1/2 ± i v^67/20). 



In the truncated case we have that eqns. (^^ and ( |90| ) 
are satisfied close to the FLRW equilibrium point, but 
now with 



1 
2^ 



570 - 1 ± V(57o - 1) + 470 (I270/5 - 1) 



(92) 

We first note that for 70 > 1/5, at least one of the 
a'' has a positive real part which leads to a change of 
stability (indeed, for 1/5 < 70 < 5/12, this equilibrium 
point is a source !). Physically, this means that in the 
truncated case S and Q in equation ( po[ ) have a growing 
mode, and hence their magnitudes increase leading to a 
breakdown in the physical model (and the time period 
for which the assumptions are valid). This breakdown is 
seen in the numerics (figures 7). As a comparison, for 
7o = 1/2, a± ~ |(1 ± 1.1); the growing mode which 
evolves approximately as e^"^, leads to a rapid increase 
in the magnitudes of S and Q and the models fail after 
a relatively short cosmological time (see figures 7) . 



VIII. NUMERICAL INTEGRATION 

From the dynamical analysis carried out in section 
there are no sources (at finite values) in the physical 
regime; that is, during the time period for which the 
various physical assumptions used here are valid there 
are no past attractors. It is reasonable to assume as con- 
ditions prevailing at the beginning of the regime we are 
interested in, that the universe is approximately isotropic 
and spatially homogeneous (e. g., almost FLRW with 
\Q\ <C 1 and E <C 1) and that the radiation com- 
ponent is dominant. These assumptions are consistent 
with current observations, and we are also assuming that 
some mechanism (such as, for example, infiation) at early 
times has driven the universe toward this configuration. 
Thus the model is close to the particular solution men- 
tioned in (ii) with (ri(r), ri(m), S, Q) given by (1,0,0,0). 
As noted above, this equilibrium point is a saddle, how- 
ever it is a 'stronger' attractor than other saddles in 
that it has more positive eigenvalues. Also, this equi- 
librium point lies in the invariant subspace X = Oj 
associated with the zero curvature Bianchi I sub-case. 
Since (il(r), f2(,„), E, Q) =(1,0,0,0) is an equilibrium 
point lying in an invariant set, if a model is driven toward 
this point in phase space it can stay an arbitrarily long 
time close to this point (i. e., the universe can spend an 
extended period close to this FLRW model). The uni- 
verse will then eventually begin to evolve away from this 
configuration and from the invariant set as time evolves. 
Therefore, we shall assume initial conditions for our nu- 
merical integration based on the fact that the universe 
starts evolving from situations close to this equilibrium 
point. 
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Initial conditions. 



Since we must have (for physical reasons) Q < and 
QE < 0, initial conditions close to ^(m)i S> Q) = 

(1, 0, 0, 0) can be given by: 



n(,)(0) - 1-e, 
S(0) = > 0, 



^™)(0) ^ e^6 > 0, 
Q(0) = < 0, (93) 



where e, S, Si, Qi are real small constants. The initial 



values of "total Omega" : fl(^tot) 
given by eqn. (|6l|) are: 



n 



(m) 



n 



(r) 



and X 



^^(tot)(0) 



1 



X(0) = <5-(S,)2, (94) 



so that the value of 6 reflects the deviation of ^(tot) 
from unity, while the deviation from the invariant set 
X = depends on S and (S^)^. 

• The factor 7. 

As discussed earlier, during the interactive range 
we are interested in, the various timescales must 
satisfy 



^ ^rcl ^ tc ^ , 



(95) 



where t„ =3/6, t„i follows from eqn . ([73|) and 
tc, t~f are given by eqns. ( |36| ) and (pTj). Consid- 
ering that (8i/3)^ « (l/3)«: a and inserting 
the constants appearing in eqns. ( P6| ) and (^, we 
arrive at the following initial values: 



logio(*ff)l* ~ 8.8, logio(ic)|» 



7.6, 



loglo(^7)l« 

thus suggesting the range 



3.4, 



0.06 < 7, < 1 



(96) 



(97) 



for the initial value 7(0) — 7^ pg[ |. 

If assuming that 7 = 70 for all times, then the 
ncar-Eckart regime can be associated with 70 < 
0(10^^), while a transient regime follows by set- 
ting 70 « 0(1) < 1. 

The constants and Qi. 



From evaluation of eqn. (53) at r = 0, the con- 
stant is proportional to the initial deviation of 
the photon entropy from its equilibrium value: 



15 Ao 



8 



(e) 



1 



^(e)(0)' 



(98) 



where we have used eqns. (11) and (||}_ so that 
Ao/<S(e) ~ 0(1), with Ao given by eqn. (^31). Both 
numbers Si, Qi are initial ratios of off-equilibrium 
and anisotropic variables (n^fc, <Jab) with respect 
to equilibrium and isotropic variables (p^'^\ 0/3). 
Since we are choosing initial conditions close to a 



saddle point associated with near FLRW conditions 
and we must assume near thermal equilibrium, then 
Si, \Qi\ must necessarily be small numbers (^ 1). 
A maximal bound on \Qi\ and |Si| can be fixed 
from CMB observations [B3| making it reasonable 
to take 



|Q,| < 0.01 -0.1, |Si| < 0.001 - 0.01, 



B. 



(99) 



However, we will comment further ahead on the 
sensitivity of the functions to these initial values. 

The constants e and 5. 

The values for e and 5 are restricted by the ratio 
of photons to WIMP's. Considering the neutralino 
as the WIMP particle with m 100 GeV, using 
the ideal gas law and eqn. (p^ yields: 



5 mc^ 



?>kgTi 



1-e f^(,)(0) 
leading to the following constraint 
5 « 1.013 e - 0.013, 



0.013, 



(100) 



(101) 



that must be satisfied by all initial conditions com- 
patible with Ti Ki 10^ K and with the observational 
constraints Vl^ ^ 0.3 ± 0.1 and h ~ 0.7. Further 
restrictions on e and 5 follow by demanding that 
r2(r) decreases and increases at the initial 

time r = 0. From the expressions for and 
O'^^^ in the differential equations, these conditions 
imply: 

-1-4S,2 < -e-6 < -2 S,|Q,;| -4 S,2, 

so that e+S > must hold, leading to the following 
minimal values of e and S : 



e > 0.0064, 6 > -0.0065. 



(102) 



The condition that fl 



(r) 



decreases at t — 0, 



together with Qi < and Si > 0, imply: 



1 



\Q^ 



< S, < -^Q^ + 4(e + 5) - (103) 



A physically plausible evolution and the range 
of validity of the models. 



It is important to specify a criterion in order to dis- 
tinguish a physically plausible evolution for the models. 
We define such an evolution by the following conditions 
that must hold all along the range of validity discussed 
previously: 



1. r2(m) increases while decreases. The transi- 

tion from a radiation- to a matter-dominated epoch 
occurs within the radiative era. However, the ra- 
tio ri(„)/f2(r-) must remain finite in all the validity 
range. 
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2. S must be an increasing and convex function, 
tending asymptotically to the equilibrium photon 
entropy given by eqn. (|35|). 

3. S'(e) oc — niust be very small. Idealy, we 
should have S'(e) > 0, though this condition might 
fail to hold as long as 5* > holds (sec for 
examples) . 



Initially, we have eqns. 
then at later stages tc 



m§, (g6|), and (g7|), but 
(the Compton scattering 



timescale) is no longer relevant, while tj and t,^i 
should overtake t^, so that the baryon-photon 
decoupling is defined as — and should occur 
at r = r„ such that T{t„) = ~ 4 x 10^ K. 

We will not be concerned with the evolution of the 
models after the radiative era, since the assumptions re- 
garding a hydrodynamic description of the radiative fluid 
break down. After the baryon-photon decoupling, an ap- 
propriate treatment of cosmic matter requires a different 
theoretical framework based on Kinetic Theory ||3^ . 

Consider a transient regime (7 « 70 « 0{1) < 1), 
together with "test" initial conditions given by eqns. ( |93| ) 
with (5 = 0, satisfying eqns. @, ( p^l ), ( p^ , and 
(|l03|), hence e = 0.0128. For the time being, will be 
taken to be two orders of magnitude smaller than Qi. 
This yields the following initial conditions lying very near 
X(0) = 0: 



17(,„)(0) = 0.0128, r!(,,)(0) = 
S(0) = 0.001, Q(0) = 
so that: n^tot){0) = 1, X(0) = 



0.9872, 
-0.1, 

-10"^ (104) 



In order to illustrate how the different variables should 
behave in a physical evolution taking place in the ap- 
propriate timescale, we integrate the system of go vern- 
ing equations (||), (H), @, and @) for ( |lo| ) and 
7 = 7o = 0.7. The results are displayed in figures 1 
and 2. As shown by figures la, lb, Ic and Id, the func- 
tions f^(r), Q, S(e) and S comply with 
conditions 1, 2 and 3 of the physical evolution mentioned 
above (we shall discuss condition 4 in section VIII-E). 
Figure 2a depicts joint logarithmic plots of the various 
timescales t-y, tc, t„, and physical time t, respectively 
given by eqns. (p6^, (p7[), ( |66| ) and (|67|), while the ra- 
diation temperature T is displayed in figure 2b. By 
comparing figures 1 and 2, it is evident that the range 
of validity of the models is roughly < r < 6, corre- 
sponding to 10^ K > T > 10^ K, with t w 10^ years 
(the physical time for the radiative era), while the tran- 
sition from radiation to matter dominance taking place 
at about r = 4 (T « 10^ K). Numerical integration 
of the governing equations for initial conditions different 
from (104) might yield important qualitative changes in 
the state variables plotted in figures 1, like r2(m), f^(r) 
or S, but not of those plotted in figures 2, such as T 
or the timescales (|6|), (^, (^), r,^i or physical time 
t. Therefore, figures 2 provide a general estimation of 
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(a) 




O.B 
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FIG. 1: A physically plausible evolution. These figures illus- 
trate the fulfillment of conditions 1, 2 and 3 for a physically 
plausible evolution given in section VIII-B. The function S(e) 
in (c) is given in cgs units and is almost equal to the equi- 
librium photon entropy. These figures (as well as those of 
figure 2) were obtained using initial conditions (102), with 
Trci = 0.7. Notice in figure lb that becomes negative 

(though small) near the initial time. This behavior does not 
denote an unphysical situation since S > holds throughout 
the evolution. 



(a) 



T = 10^ K 



T = 4 X 10^ K 




FIG. 2: Range of validity of the models. The range < r < 
6 corresponds to the temperature range of the radiative era, 
between Ti — 10^ K and the baryon-radiation decoupling 
temperature 4 x 10^ K (roughly at r = 5.5). A logarithmic 
plot of the various timescales used in the paper (in seconds) 
are depicted in figure 2b: the Thomson and Compton scat- 
terings {t.y, tc), the Hubble time t^j and physical time t. 
Notice how the decoupling temperature occurs at the same 
r as the condition t^ = t„, while radiation-matter equal- 



ity, n 
(roughly 



('■) 



^(m) (see figure la) corresponds to T - 
= 4.6.) 



10* K 



the range of validity of the models for a wide range of 
initial conditions. We will consider more general initial 
conditions in the following subsection. 
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C. Sensitivity to deviations from x — 0. 

We will test the effect of initial deviations from the in- 
variant set X = on state functions obtained by numer- 
ical integration of the governing equations. We consider 
initial conditions as in (93), keeping S^, and Qi fixed, 
but now we take 5 ^ Q: 



rj(„)(0) = 0.0128-5, 
E(0) = 0.001, 
%ot)(0) = l-<5, 



r2(^)(0) ^ 0.9872, 



Q(o) 
x(o) 



-0.1, 

5 - IQ- 



(105) 




Since 5 can be given in terms of e by eqn. (101) and 
e = 0.0128 corresponds to (5 = 0, testing values of e 
near e = 0.0128 determines the initial deviation from 



FIG. 3: Sensitivity to x(0) = 0, positive curvature. Figures 



display Q, 



and S/S(e) for 



(tot) 



1 and X = 0. Notice that 5 can be positive initial conditions (103) with —0.01 < 5 <Q\ numerical values 



or negative, respectively, for e > 0.0128 or e < 0.0128, 
so that x(0) = if (5 = Sf = 10-^ while x(0) is 
positive/negative ii 5 > = lO"*^ or (5 < = 10"^ 
(though, from eqn. (|6l|), curvature has opposite sign to 
x). We integrate nu meric ally the governing equations 
for initial conditions (105), 7 = 70 = 0.7 (transient 
regime) and for various values of 5. The resulting forms 
of r2(m), f2(r), f^(tot) I and S are respectively plotted in 
figures 3a to 3d (for positive curvature X < ) and in 
figures 4a to 4d (for negative curvature x > 0)- These 
figures clearly show that a physical evolution is only pos- 
sible for initial conditions that deviate very slightly from 
x(0) = (less than w 10~^-^), leading to orbits that 
remain very close to the invariant set x = 0. If we fix 
e and 5 (for any combination of values compatible with 



( |101| ), ( |102| ) and ( |103| )) and vary S^, so that deviation 
from x(0) — is governed by , we obtain exactly the 
same behavior displayed by figures 3 and 4, leading to the 
same conclusion: a physical evolution is only possible for 
initial conditions for which |x(0)| % 10^^'^. 



of 5 appear next to each curve. Q. 



(m) ; 



(r), 



and Q, 



(tot) 



branch upward and S/S(e) downward; the branching time r 
decreases with increasing \S\. The branching up corresponds 
to 0^0, marking the re-collapsing stage of the models. 
A comparison of figures (a), (b), and (c) with (d), and with 
figures 2a and 2b, shows that a physically plausible evolution 
for the duration of the radiative era (0 < r < 6) and for 
the appropriate temperature ranges is only possible if \5\ < 





D. The near-Eckart regime 

We examine the near-Eckart regime by assuming 7 — 
7o = 0.001 <C 1, together with initial conditions ( 104 ). 
As shown in figure 5a, the functions Q{t) and S{t) 
clearly have the form (|76|), indicating a quick relax- 
ation in terms of an abrupt exponential decay (in about 
r « 0.01). Figure 5b shows how the equilibrium entropy 
S'(e) tends to a constant value in about t « 0.02, a 
longer time than the relaxation of S/S^^y^ thus indi- 
cating that the effective relaxation time is provided by 
S'(e), and 5^^-! > holds for all the evolution, in agree- 
men t with the "local equilibrium hypothesis" (see sec- 
tion [yiA] ). Figure (5c) reveals how (5/8) Q/S 70 = 
0.001, so that Eckart's transport equation {P + 2ija = 
with 77 given by eqn. (p9[)) is approximately valid once 
the quick relaxation is over. A good approximation to 
the relaxation time in a near-Eckart regime is given by 
the assumption (^ with 7 = —(5/8) 70 (Q/S) and 
7o = 0.001. Figure 5d illustrates that this is a correct 



FIG. 4: Sensitivity to x(0) = 0, negative curvature. Ana- 
logues of figures 3 for 5 > 0. A similar branching of the func- 
tions fi(r), ^(tot), 'S'/5'(e), is obscrved cveu though the 
models do not re-collapse. Again, the function S'/S'fe) sat- 
isfies the conditions for a physically plausible evolution only 
if initial conditions are given by 5 < 10~*'^. Notice that the 
^(m){T) curves in (a) with S ~ lO^'''^, bend downward, so 
that for r ^ 6 they tend to the currently accepted values of 



present CDM abundance: Q 



(m) 



0.3. 



assumption, since the obtained r,ei overtakes in the 
very short period that coincides with the duration of the 
relaxation process (r « 0.02). 

The plots of the functions ri(-„j)(r) and fl(^j,-^(T) are 
identical to those that would have resulted i n th e tran- 
sient regime had we chosen initial conditions ( |l04| ) and a 
much larger value of 70. The function E(r) is affected, 
becoming almost constant for very small 70. This is 
reasonable, since eqn. (pq) does not contain Q, while 
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(a) 



(bV 



the numerical valu(3 df 



those shown in figy ir( ps 1 to d, being l^h^us unaffected by 

not Vpry realistic, we test now 
that follows from equations ( |8^ ) 
lIH. Ignoring the curvature term, 



Since 7 constcini 
the expression for 
and ( p^ ) in section 
-^w e hav e that fM)/"/ 



■+Mo(l + Mo) 



Mo 



18, 



(40§)- 



FIG. 5: "Near-Eckart" regime. Functions S/S[e), S(e)/kg, 
and — (5/8)((5/E) correspond to initial conditions (104) with 
70 — 0.001. S/S(_f,) in (a) relaxes much faster (r « 0.002) 
than in figures Id, 3d or 4d, associated with the transient 
regime. The decay of 5(e) in (b) takes longer (r ~ 0.01) 
than that of S, hence 5(e) is an adequate entropy function 
in agreement with the "local equilibrium" hypothesis charac- 
teristic of the Eckart regime (5(e) > holds throughout the 
evolution). Panel (c) shows how -(5/8)(Q/E) ^ 70 = 0.001, 
indicating that Eckart 's transport equation is approximately 
valid for r > 0.02. In panel (d) we used initial conditions 
(104) and 7 = -(5/8)(Q/E)7o with 70 = 0.001 (instead of 
7 = 70 = 0.7), which yields an excellent approximation to the 
relaxation parameter of a near-Eckart regime, overtaking 
at the relaxation timescale t « 0.01; this timescale however, 
is too short in comparison with the relaxation timescale of the 
Compton scattering also shown in (d), hence this and other 
radiative processes must be studied within a transient regime. 



eqns. ( |58D and (q9|) do contain this function, but \Q{t)\ 
and Q' decay very fast becoming almost zero for most 
of the time range and so the differential equations for the 
functions f2(m) and Q(^r) (but not S) are practically 
unaffected. The effect of varying 70 (as long as we have 
small values, < 0.01) is simply to make the decay of 
Q and S'f^^^ slightly more or less abrupt (depending on 
whether 70 is smaller or larger than 0.001) and has no 
noticeable effect on numerical curves of other functions. 



It is evident that choosing a sufficiently small /iq yields 
7i ~ 5/io/8 ^ 1, an initial value characteristic of a 
near-Eckart regime. As expected, numerical tests with 
/iQ ^ 1 lead to practically identical curves as those cor- 
responding to the near-Eckart regime (figures 5a, 5b and 
5c). However, for /iq > 1/2; the factor multiplying 



ri(„j) in the denominator of ( |106| ) can become negative 
for some t, thus opening the possibility that the de- 
nominator might become small and so 7 might increase 
to « 0(1) as T reaches later times related to the 
baryon-photon decoupling era. In order to explore this 
possibility, we in tegrated the dynamical system for initial 
conditions (|104| ) and under the assumption of 7 given 
by eqn. (|106| ), with various values of /iq > 1/2. As 
shown by figure 6a, the choice /io > 15 leads to 7 
diverging as (approximately) r ^ 7, this leads to t„i 
overtaking (figure 6b) around t « 6.8. It would 
have been nicer to have t,„i overtaking at an ear- 
lier time (say r « 5 — 5.5) as required by condition 4 
for a physically plausible evolution, but we feel that the 
form of r„i associated with (106) is a reasonable ap- 



proximation to a physical relaxation time that acts as an 
"effective" relaxation time for the radiative era. Finally, 
another consequence of dealing with a more reasonable 
form for t„i is the fact that ^(e) oc — E Q being posi- 
tive and very small (condition 3 of a physically plausible 
evolution) is better satisfied than in the case of constant 
7o (compare figures lb and 6c). 



E. Testing the relaxation times numerically. 



The assumption 7 = 70 < 1 leads to a reasonable re- 
laxation time only in the earlier stages of the evolution, 
while a choice 7 = 70 > 1 might work for later condi- 
tions (near matter-radiation decoupling) but not for ear- 
lier times. This can be appreciated in figure 2b, since the 
relaxation times obtained for different values of 7 = 70 
would be curves parallel to . As we mentioned before, 
7o controls the rate of decay of Q and S, but it is 
still interesting to check if other functions are sensitive to 
changes in the numer ical constant value of 79. Assum- 
ing initial conditions (104) and integrating the dynamical 
system for various values of 70 > 0.1 (thus excluding 
the non-transient zone) shows that the other functions 
(such as ri(m): ^(r) or E) are essentially identical to 



F. The truncated equation. 



Considering initial conditions (104) and 7 = To — 0.7, 
the integration of the governing equations (pq), (|5S 



I), and (|86[) yield the curves for Q, E, and 



ri(r) depicted in figures 7a, 7b and 7c. These figures re- 
veal that a physical evolution fails to occur, as the grow- 
ing modes predicted by the qualitative analysis clearly 
emerge, making all these functions undertake an unphys- 
ical growth. We te sted other values of 70 and 7 
given by eqn. (106), as well as other initial conditions 
and obtained very similar curves to those of figures 7a- 
c, all failing to comply with the criteria for a physically 
plausible evolution. 
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FIG. 6: Dynamical relaxation time. Panel (a) displays 7 ob- 
tained from integrating the governing equations for the form 
given by (106), the numbers next to each curve correspond to 
the chosen numerical values of A^o; 7 diverges for fio > 15. 
Panel (b) displays the corresponding form for t^^i which 
overtakes t^^ for r ^ 6.8. Panel (c) shows S(^e)/kg for 
fio = 200. A comparison with figure (lb), shows that requir- 
ing S[s)/kg to be small is better satisfied with 7 given by 
eqn. (104) than with a constant value 7 = 0.7 (as in figure 
lb). 



G. A baryonic scenario without CDM 

If we had considered only the radiative fluid (baryons, 
photons and electrons) as the matter source of the mod- 
els, then the bulk of the rest mass energy density would 



FIG. 7: The truncated transport equation. Curves obtained 
by integrating the governing equations with initial conditions 
(104) and 70 — 0.7, but using eqn. (83) ("truncated" trans- 
port equation) instead of eqn. (64) (full transport equation). 
The exponential growth terms mentioned in section VII lead 
to an unphysical evolution characterized by negative 5(e) 
and S during all the evolution range. 



ri(r) larger, since baryons have less rest mass density 
(by one or two orders of magnitude) than WIMP's and 
so it should take longer for baryons to dominate over 
radiation. The numerical curves that result are as ex- 
pected intuitively, with ri(„i), f^(tot), and S/S^^-f hav- 
ing very similar forms as the curves of figures 3 and 4, 
with ri(r) decreasing slightly slower than in the case 



(m) 



(™) would have 



(m) 



have been due to the baryons, so that n 
been identified with rib instead of n„ and Q 
would correspond to the baryonic f2b- It is interesting 
to examine numerically the consequences of this "bary- 
onic scenario" . Considering the baryon-photon number 
ratio i n eqn . (|l5|) , the baryonic scenario implies replacing 
eqns. (|10C|) and (101) by 



-S 



1 



0.002, 



6 = 1.002 e- 0.002, 



(107) 

while initial conditions are then given by eqns. ( 104) and 
( |105| ), but S — now corresponds to e = 0.00199 in- 
stead of e = 0.0128. Intuitively, we do not expect a 
mayor qualitative change in the resulting graphs, though 
it is reasonable to expect that fl^m) will be smaller and 



with WIMP's. Since the obtained curves for Q 
the baryonic scenario are so close to those obtained in fig- 
ures 3 and 4 for the case with WIMP's, these curves yield 
baryon abundances that are clearly incompatible with the 
bounds placed by cosmic nucleosynthesis (fit ~ 10"-^). 



IX. DISCUSSION AND CONCLUSION. 

We have studied a class of dissipative Kantowski-Sachs 
models describing the cosmological evolution during the 
radiative era characterized by radiative processes involv- 
ing baryons, electrons and photons, considered as a single 
dissipative "radiative" fluid. We also assumed the pres- 
ence of CDM, in the form of a non-relativistic gas of 
WIMP's (lightest neutralinos) . Although this gas does 
not interact with the radiative fluid, it provides the bulk 
of the rest mass energy density and thus it strongly in- 
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fluences the dynamics of the models and the resulting 
values of cosmological timescales, such as t^. On the 
other hand, the radiative fluid provides the bulk of ther- 
mal and dissipative effects, related to the rate of change 
and relaxation of the radiation entropy to its equilibrium 
value. 

After defining new normalized variables ^(rn)i ^(r), 
S, Q, a set of evolution equations has been derived based 
upon appropriate thermodynamical laws and equations 
of state. The qualitative and numerical study of these 
evolution equations has clarified various aspects of the 
dynamical behavior of the models, their physical viabil- 
ity, as well as a peculiar sensitivity to certain initial con- 
ditions related to deviation from the invariant set x = 0. 
We discuss below the main features emerging from pre- 
vious sections. 

The definition of the phase space variables fl(^r) : ^^(m) : 
S, Q leads in a natural way to expressing the relaxation 
time T,„i as proportional to t^^ =3/9 (see equations 
( |68| ) and (73)). The understanding of the relaxation pro- 
cess can be accomplished by studying the effect of dif- 
ferent choices of the proportionality factor, 7(r) > 0, 
on the exponential decay of the dissipative stress (related 
to Q) and of the photon entropy 5' to its equilibrium 
value 5(2). We have identified a "near-Eckart" regime 
if this decay is abrupt (7^1, and becoming instan- 
taneous in the limit 7—^0, so that r„i — > 0), while 
a "transient regime" (7 w 0(1) < 1) can be associated 
with a slower decay. Both the near-Eckart and the tran- 
sient regimes are compatible with a physically plausible 
evolution. The difference between the two regimes is the 
timescale of their relaxation process: for the transient 
regime this timescale can be comparable with the dura- 
tion of the radiative era, for the near-Eckart regime it is 
much shorter (about eight orders of magnitude in phys- 
ical time). This is well illustrated by the differences in 
evolution timescales between figures 5 and figures 1 to 
4. Comparing figures 2b and 5d, it is evident that the 
relaxation timescale of the near-Eckart regime (« 10^ 
sec) is much shorter than that of the Compton scatter- 
ing (« 10^^ sec). Therefore the near-Eckart regime yields 
a relaxation that is too swift and so it is inadequate to 
examine the two main radiative processes of the radiative 
era: the Compton scattering and (more so) the Thomson 
scattering. It is important to mention this fact in view of 
recent claims that a transient theory of irreversible 
thermodynamics is not really necessary (see |]36| for a 
comprehensive discussion). It is evident that the cos- 
mological study of radiative processes in pre-decoupling 
times needs to be accomplished with a transient regime. 

As revealed by figures 1,2,3 and 4, a physically plausi- 
ble evolution is possible for all the duration of the radia- 
tive era for initial conditions given by an initial state very 
close to x(0) = 0, hence lying very close to the invari- 
ant set of zero curvature x = 0. All models complying 
with a physical evolution begin their evolution near the 
equilibrium point {i^(r)j ^(m)i 5], Q) = (1,0,0,0), a 
saddle with positive eigenvalues (i. e., stable) associated 



with the FLRW sub-case of the models, and the pro- 
ceeding evolution remains very close to the invariant set 
X = 0. This is an extremely interesting feature of these 
models, as it relates a geometric property of KS solu- 
tions with the constraints imposed by the physics and by 
observational evidence, since recent data from CMB ob- 
servations indicates fl 



of and fl 



(•jgj^ « 1. By looking at the curves 
tot) with 6 < 10^^ in figures 4a and 
4c (negative curvature), it is evident that for r > 6 
these curves decrease from their values r2(m) ~ 1 and 
^(tot) ^ 1 around t = 6. Had we plotted these curves 
for larger values of r, , extending to the present era 



(r w 15), we would have obtained O 



(m) 



n 



{tot) 



0.3, 



in agreement with the currently accepted value of r2(„i), 
but not of ri(fot)- Of course, the estimated contribution 
to r2(tot)j today, for non-relativistic matter (CDM plus 
baryons) is only « 0.3, with the remaining two thirds of 
the critical density probably related to a A-type "dark 
energy" interaction whose precise nature and properties 
are still uncertain. However, this discrepancy with re- 
gards to ^l(tot) today is not surprising since we did 
not consider any A-type interaction, and so it does not 
affect our results. Also, the models we are considering 
are only valid for a specific range of cosmological times: 
10'^ ^ -2^ ^ 10^: in which this "dark energy" would likely 
not have been dominant. Still, the close link between a 
physically plausible evolution and ri(tot) near unity is 
remarkable. 

It is interesting to compare our results to those re- 
ported previously |0 dealing with the perfect fluid sub- 
case of the KS models examined in this paper (though, 
these models did not consider a CDM component). 
As reported in js^ , there are numerical solutions in 
which both matter and radiation normalized densities, 
r2(„j), r2(r), decay to zero as the models re-collapse and 
approach a crunch singularity. By looking at the forms 
of ri(m) and ri(r) m figures 4, it is evident that such 
evolution is similar to that depicted by curves associated 
with initial conditions 6 > 10~^. However, the evolution 
that results from these initial conditions fails to comply 
with our physical criteria, since the entropy S is no 
longer a convex function for all of the time range (see fig- 
ure 4d) and starts decreasing at too early times. In the 
perfect fluid case, these examples satisfy an appropriate 
equation of state and all of the energy conditions and 
also the photon entropy is simply S'(e) and is constant 
for all times, hence there is no physical reason to discard 
these curves (other than remarking that such behavior of 
r2(^), fl(^r) is not observed in the real universe). How- 
ever, for the dissipative source under examination here, 
S and S(^e) are not constant and the conditions for a 
physical evolution are more stringent, hence we can ap- 
ply clear physical criteria to discard these perfect fluid 
"strange" cases. 

For initial conditions near the equilibrium point (ri(r), 
S, Q) = (1,0,0,0) (as in equation (p^), the 
numerical curves of the "equilibrium variables" ^^(m) 
and ri(r) are not affected by the choice of 7 = 70, the 
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constant proportionality factor between r^oi and . 
This is a consequence of the fact that for these initial 
conditions the values of S, Q and S — ^(e) remain 
small during all of the evolution. Therefore, the evolu- 
tion equations for fl(^„i) ^nd f2(r) (equations ( p6| ) and 
(^)) are practically unaffected by the presence of S 
and Q, and so are insensitive to the rate of transiency 
given by 7 constant. Since E and Q govern the devi- 
ation from the FLRW equilibrium point, this decoupling 
of f2(„i) and ft(^j.) from E and Q is then a consequence 
of the evolution of the system always remaining close to 
thermal equilibrium. We tested this behavior for a par- 
ticular form of 7 (eqn. (|106|), section VI B): the curves 



for i^{r), E, and Q are practically identical with 

those that follow from choosing 7 = 0.7 in figures 1 
to 4. However, as shown in figures 6a and 6b, for large 
values of /iq defined by (|8^), the obtained relaxation 
time r,„i behaves similarly to what one would expect of 
a relaxation parameter for the radiative era. Although 
it has become common practice to simply equate r„i 
with a microscopic interaction time, like or tc, the 
relaxation time is not a microscopic but a mesoscopic or 
even macroscopic quantity (though it must be qualita- 
tively analogous to interaction times [Q). Since it can 
be extremely cumbersome to evaluate t^^i , it is useful to 
have a concrete example where this relaxation parameter 
can be adequately approximated by the same dynamical 
equations associated with the models. 

Finally, by means of qualitative arguments supported 
by the numerical analysis, we have shown in section VII 



(figures 7) , that the truncated equation does not comply 
with a physically plausible evolution. This is an impor- 
tant result, since we have found a concrete example in 
which a truncated transport equation leads to unphysi- 
cal evolution of dissipative fluxes. Although this conclu- 
sions strictly applies to the KS models under considera- 
tion, we must point out that one should be very cautious 
when applying these equations to other models and other 
equations of state. 

A possible extension of this work would be to con- 
sider instead of CDM other forms of dark matter, such 
as "warm" dark matter (WDM) or axions. Another pos- 
sibility is to include, together with dark matter, a scalar 
field associated with "dark energy". Another route to 
generalize the present work is use a class of metrics asso- 
ciated with a geometry that is less restrictive than KS, 
for example the non-static spherical symmetry (perhaps 
under the assumption of self-similarity). We regard the 
present analysis of the Kantowski-Sachs models as a first 
step toward an understanding of the dynamics of cosmic 
matter in more general and physically motivated inho- 
mogeneous models. 
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